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This paper explores various interfaces between the class of generalized gamma convolution (GGC) 
random variables, Dirichlet process mean functionals and phenomena connected to the local time and 
occupation time of p-skew Bessel processes and bridges discussed in Barlow, Pitman and Yor (1989), 
Pitman and Yor (1992, 1997). First, some general calculus for GGC and Dirichlet process means func- 
tionals is developed. It then proceeds, via an investigation of positive Linnik random variables, and more 
generally random variables derived from compositions of a stable subordinator with GGC subordinators, 
to establish various distributional equivalences between these models and phenomena connected to local 
times and occupation times of what are defined as randomly skewed Bessel processes and bridges. This 
yields a host of interesting identities and explicit density formula for these models. Randomly skewed 
Bessel processes and bridges may be seen as a randomization of their p-skewed counterparts, and are 
shown to naturally arise via exponential tilting. As a special result it is shown that the occupation 
time of a p-skewed random Bessel process or (generalized) bridge is equivalent in distribution to the 
occupation time of a non-trivial randomly skewed process. 

1 Introduction 

This paper explores various interfaces between the class of generalized gamma convolution (GGC) 
random variables, Dirichlet process mean functionals and phenomena connected to the local time and 
occupation time of skew Bessel processes and bridges discussed in Barlow, Pitman and Yor (1989), 
Pitman and Yor (1992, 1997b). The latter concepts are also related to the now ubiquitous two- 
parameter Poisson-Dirichlet models, see for instance Pitman and Yor (1997a) and Pitman (2006). 
Specifically, we take the viewpoint as in James (2006), that GGC random variables can be represented 
as linear functionals of Gamma processes and hence there is a strong link to the distributional 
theory of Dirichlet mean functionals developed in Cifarelli and Regazzini (1990). Here, we first 
develop two types of distributional results, one is a simple but quite useful result involving scaling 
general Dirichlet mean functionals by beta random variables, and the other is a form of gamma tilting 
calculus. We then investigate the class of positive Linnik random variables(see Devroye (1996, 1990)), 
establishing new distributional identities, and a link to local and occupation times of skew Bessel 
processes and bridges. As a nautral extension, we investigate compositions of a stable subordinator 
and a subordinator which has a GGC law. We use this to develop additionally a scaling calculus 
for (generalized) Dirichlet mean functionals. These points then lead to the development of a series 
of interesting results, including explicit density formula, for what we call randomly skewed Bessel 
processes and bridges. As a special case, we establish equivalences in law between occupation times 
of p-skewed Bessel processes and bridges and non-trivial randomly skewed processes. Some recent 
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work related to our line of investigation include, James (2006), James, Lijoi and Priinster (2006), 
Bertoin, Fujita, Roynette and Yor (2006) and Fujita and Yor (2006). With some further reference 
to these latter two works, we note that a key result in this manuscript is to show that the random 
variable X a = S a /S' a , where S a and S' a are iid a-stable random variables is indeed a Dirichlet mean 
functional. The explicit distribution of X a was obtained in Lamperti (1958) and, as shall be shown, 
this fact combined with the correspondence to Dirichlet mean functionals allows one to easily obtain 
explicit densities for a host of models connected to Bessel phenomena. For some other references 
relevant to skew Bessel processes see Bertoin and Yor (1996), Kasahara and Watanabe (2005) and 
Watanabe (1995). 



2 The relationship between GGC, Gamma processes and 
Dirichlet mean functionals 

A positive infinitely divisible random variable Zg is said to be a generalized Gamma convolu- 
tion (GGC) if for positive A one can write its Laplace transform in the form 

E[e- AZ »] = e~ M W 

where for a sigma- finite measure v on (0, oo), 

0(A) = / log(l + Xx)v{dx) = / / {l-c-^s^e-'^vidx). 
Jo Jo Jo 

v is chosen so that -d(X) < oo, and is often referred to as a Thorin measure. For precise conditions 
on v sec Bondesson (1992). The corresponding Levy measure of Zg is, 

/>00 

6s- 1 \ e~ s ' x v{dx). 
Jo 

We will say Zg is GGC(#, v). Naturally every GGC can be used to generate a subordinator which 
we denote as Zg(t) for each t > 0. We will also be interested in the special case where v = H 
is a probability measure, in that case we write GGC(9,H) and will sometimes refer to Zg as an 
FGGC or finite GGC, a term we used in James (2006). As in that work, we will exploit a close 
connection between GGC, and Gamma processes, and in the case of FGGC, its connection with the 
Dirichlet process. [See Ferguson (1973) for the Dirichlet process and its use in Bayesian nonparametric 
statistics]. Now let Tg v denote a gamma process on (0, oo) with shape parameter (0v), this is a 
completely random measure whose law is characterized by its Laplace functional for some positive 
function g as 

E [ e -ArV(g)j = e -e/ °»log(l+Ag(x)) ! ,(dx) ! 

where Tg u {g) := J °° g{x)Tg l/ (dx). Hence setting g(x) = x, we see that every GGC(#, v) random 
variable can be written as a mean (linear) functional of Tg v . Let (Jk,e) denote the ranked jump 
sizes of a subordinator based on a Gamma(#) distribution, such that X^fcLi Jk,e — Gg. Now recall 
that for each 9 > 0, Gg, which is a Gamma(6>) random variable, is independent of the sequence 
of ranked probabilities (Pk,e = Jkj/Gg). The law of the sequence {Pk,e) is sometimes referred to 
as the Poisson-Dirichlet distribution and, using notation from Pitman and Yor (1997a), is denoted 
PD(0, 0). These points imply that 

xTg v (dx) = GgMg{v) 

where Mg{v) = J2T=i Pk,eWk, where the (Wk) are independent of (Jk,e) an d are the points of a Pois- 
son random measure with mean intensity v. Furthermore, the distribution of Mg(y) is characterized 
by its Cauchy-Stieltjes transform of order 8, which is 

E[(l + \Mg(v)y e ] = E[e- XGeMsiv) ] = c-^ (A) . 
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Hence, importantly, it follows that Zg = GgMg(v). Now when v — H we may define a Dirichlet 
process with shape parameter 9H as Po,g{-) — Y^k=i Pk,e5z k (•)> where the (Zk) are iid with common 
distribution H. Then a Dirichlet mean functional is defined as 

r-OO 

Mg{H) = / xP 0i g{dx). 

Jo 

Hence, if Zg is a GGC(0, H) then Zg = GgMg(H). 

Remark 1. Note of course that when v has infinite mass the use of 9 is somewhat redundant. 
However, it will play an important role as the decomposition GgMg(y) always makes sense. 

Remark 2. In Bondesson (1992) the Levy exponent of a GGC is represented in a slightly 
different, but certainly equivalent, way as 

/>oo 

0(A) = / log(l + \/y)U{dy). 
Jo 

where U is some sigma-finite measure. In this work, we are using a representation which more closely 
matches that used in the literature on Dirichlet process mean functionals. 

Remark 3. See Vershik, Yor and Tsilevich (2001) and James (2005) for some uses of the 
decomposition GgMg(H) 

Remark 4. Throughout, for < a < 1, we will define S a to be a unilateral (a)— stable random 
variable having the Laplace transform 

E[ e - AS «] = c - xa . 

This random variable is well known to be a GGC with an infinite Thorin measure. See example 3.2.1 
of Bondesson (1992). 

3 Some calculus for GGC/Dirichlet mean functionals 

Suppose that X has distribution H, and define the function 

$(t) = / logflt - x\)I{t + x)H(dx) = Epogfli - X\)I(t + X)} 
Jo 

furthermore, letting H(t) — J Q H(dx), define 

A e (t\H) = -sm(TT9H(t))c^ t \ 

Cifarelli and Regazzini (1990, 1994), see also Cifarclli and Mcllili (2000), apply inversion formula to 
obtain the distributional formula for Mg(H) as follows. For all 9 > 0, the cdf can be expressed as 

(1) C {x-t) e -^ e {t\H)dt 

Jo 

provided that 9H possesses no jumps of size greater than or equal to one. If we let /m 8 ('\H) denote 
the density of Mg(H), it takes it simplest form for 9 = 1, which is 

(2) fM 1 (x\H) - A^xlH) = -sm(nH(x))e-^ x \ 
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Formula for 9 > 1 are also available from Cifarelli and Regazzini (1990,1994) and Cifarelli and 
Mcllili (2000). However, we shall not explicitly use this and instead give an expression for the 
density which holds for all 9 > 0, recently obtained by James, Lijoi and Prunster (2006), as follows, 

(3) f Me (x\H)= f (x-tf-'deit^dt 

Jo 

where 

d e {t\H) = j t ^n{n9H(t))c- e ^ l \ 

We do point out, that except for some of the very recent work we mention here, there are very 
few examples of mean functionals where the explicit density has been calculated. This is due in 
part to the fact that it is not necessarily obvious how to write Ag(t\H) in a nice form, or otherwise 
how to calculate $(x). Our results in this section provide some additional tools to obtain more 
explicit expressions and also point to apparently unknown inter-relationships between different mean 
functionals which occurs through a tilting operation. 

The first result, in Theorem 3.1 below, shows that multiplying a mean functional by a certain 
beta random variable may lead to more simplified expressions. Additionally, we point out that this 
result was constructed in an effort to better understand the result in Proposition 15 of Pitman and 
Yor (1997b) and, as we show, paves the way for other interesting results related to occupation times 
and local times of Bessel processes. 

Remark 5. First, we introduce a bit more notation. The notation G a , for a > 0, will denote 
a Gamma random variable with shape a and scale 1, with law denoted as GAMMA(a). Throughout 
we will use the notation B a ^ to denote a Beta random variable with parameters (a, b), and write 
BETA(a, b) to denote this law. Additionally, unless otherwise specified, if for random variables X 
and Y we write the product XY, it will be assumed that X and Y are independent. Additionally 
the notation X' means that X' is an independent random variable with the same distribution as X 

Theorem 3.1 Let X denote a random variable with distribution H. For < p < 1, let Y p de- 
note a Bernoulli(p) random variable independent of X. Then the distribution of XY p is given 
by H^(dv) := pH(dv) + (1 - p)S (dv). Now define $(t) := E[log(|t - X\)l(X ^ t)}. Hence 
$fr)(t) := Epogflt - XY p \)l(XYp ^ t)} = p<$>{t) + (1 - p) Iog(t)I(t + 0). Then for 9 > 0, 



(4) M e (H^) £ M pe (H)B w . 



■p) 



and G e M e {H^) = Gg p Mg p (H). Hence GGC(6, H^)=GGC{9p, H) for all 9 > Oand < p < 1. As 
special cases, 

(i) Mi(H^) = M p (H)B Ptl - p , with density, for < p < 1, 

(5) -aT (1 - p) sin(7rp[l - H{x)])e-^ {x) 

TT 

(ii) Mi(ffW) = M 1 (H)B 1 (i-p) . 



Proof. First note the distributional identity, Gg p = GbBq p qci-p), for all 9 > 0. Then for (4) 

it suffices to check the equivalence GeMg(H^) = Ge P Mg p {H). But this is obvious upon taking 
Laplace transforms. The density in (5) follows as a special case of (2). Note additionally the identity 
sm(-K[pH{x) + l-p}) = sin(7r[p(l - H(x)]), for < p < !.□ 
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3.0.1 Remark about obtaining the distribution of a FGGC subordinator for all time 
points 

It is well known that if {Z(t),t > 0} is some subordinator, even if we have an explicit expression for 
its density or distribution at time say 1, i.e. the random variable Z(l), it is not obvious how to find 
explicit densities or distribution functions for general Z(t), at fixed time points t. Of course, by the 
convolution properties of infinitely divisible random variables one can describe the law of a random 
variable for Z(t) for t > 1, in terms of a linear combination of independent random variables based 
on {Z(s) : s < 1}. Hence, from this point of view, one may simply concentrate on obtaining explicit 
expressions for the density of random variables for Z(t), where t is some fixed value in (0, 1]. In the 
case where {Zg(t);t > 0} is a FGGC subordinator, one may obtain an explicit description of its 
density for each fixed < t < 1, provided that $ has a tractable form, using Theorem 3.1. Formally 
if < 6t < 1, it follows that 

Z e (t) = Ge t M et {H) = GiB w>1 _« t M 9t (ff) = G\M\{H^) 
where Mi(H^) has a density (of non-integral form), 

-x 81 - 1 sin(7T0t[l - H{x)])e~ et ^ x \ 

The point here is that density of Mi(H^) generally has a simpler description than that of Mg t (H). 
Naturally, these models agree when Ot = 1. 

Remark 6. These points and Theorem 3.1 are connected to an expression for a GGC(/3, H) 
density for < /3 < 1 given in Bondesson (1992, p. 37-38). Relationships to beta random variables 
are not noted there. 

3.1 Some basic exponential tilting calculus for Gamma mixtures 

Here we develop distributional results for exponential tilting of scale mixtures of Gamma random 
variables. We then specify this to the case of exponential tilting of GGC models which always can 
be written in this form. An important point is that we are able to establish a clear distributional 
link between a Mg(v) and a corresponding random variable obtained from tilting GgMg(v), which, 
as a by-product, has new implications for the study of Dirichlet mean functionals These results will 
play a fundamental role throughout the text. 

Theorem 3.2 Suppose that W = GgM where M and Gg are independent positive random variables 
and Gg is Gamma(#). Then for b and c positive numbers, the random variable, Wg with density 
e -c/ bw b - 1 f w (w/b)/E[e- cG <> M ] satisfies, 

Wg = {b/c)GgYg, c 

where, 

(i) Yg tC is a random variable with density 

fy(y)[(i-y)} e 

E[(l + cMY ] ' 

where Y = cM/(cM + 1), and E[(l + cMy 9 ] = E[ e - cG « M ]. 

(ii) Equivalently the distribution ofYg tC is the distribution ofcM/(cM+ 1) taken with respect to 
the density (cm + l)~ e f M (m)/E[ e - cGsM ]. 
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(Hi) Conversely f M {x) = (1 + xf 2 f Y$1 ( + M) 9 ] 

Proof. Let us proceed by checking Laplace and Cauchy transforms. By direct argument, work- 
ing with the density the Laplace transform of Wg evaluated at A is, 

E[(l + (c + Xb)M)- e ] 



E[(l + cMY°] 

which, for c 7^ 0, may be re- written as 

E[(l + c(l + \(b/c))My 6 } 



E[(l + cMY u ] 

Now for simplicity we set b/c = 1 and note that the Laplace transform of GgYg jC is, by definition of 

Y e ,c, 

cx \~ (1 + cx)~ f M (x) ^ 



E [U + AY„,r«] = f (l + A-f^ 



E[(l + cM)" 



Simple algebra completes the result. □ 

The next two results are specialized to GGC. 



Proposition 3.1 Let L be a GGC{9 1 v) random variable with density denoted fh- Hence L = 
GgMgiv). Now define a tilted random variable L having density proportional to e^ c / hw b^ 1 fw(w /b) 

fore > andb > 0. Denote this law of L as GG& b '°}(6,v). It then follows that for c > L = (b/c)L c 
where L c is GGC^ c ' c \0,u). 

(i) The Levy measure of L is 

1*00 

s -l e -(c/6)» / e -«/»r I/ (d r ) 







(ii) When c = 0, L = bL. Otherwise for c> 0, L = {b/ c)GgMg(v^). The Thonn measure 

is determined by v as indicated by the following expression for the Levy density of GgMg{v^ c ^), 



POO POO 

s- 1 e- s ^v{dr) = s- 1 e- s,y v {c - c \dy) 
Jo Jo 

(in) Referring to Theorem 3.2, M e (u^) = Yg^ c , for M = M e {v). 



PROOF. The first statement is just a standard result for exponentially tilting infinitely divisible 
random variables coupled with the specific form of a GGC Levy measure. Statement [(ii)] is just 
algebra, which yields importantly the representation GgMg{v^ c ^). Statement [(iii)] then follows 
from Theorem 3. 2D 

We close this section with an important variation in the case where v is a probability measure. 

Proposition 3.2 Let L be a GGC{9 1 H), then there exists a random variable X with distribution 
H and a random variable A* = cX/cX + 1 on [0, 1] such that if L c is now it is 

equivalently GGC(9,Q C ), where Q c is the distribution of A* derived from H and is a special case of 
j,(cx) j n ^ er wor d s th e l,£ V y measure of L c may be written as 

Additionally the following distributional relationships hold. 
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(i) Suppose that the density of Mg(H), say fM e (:\H) is known. Then the density of Mg(Q c ) is 
expressible as 

(ii) Conversely, if the density of Mg(Q c ), fM e (-\Qc), is known then the density of Mg(H) is given 
by 

f Me (x\H) = (1 + x) e - 2 f Me ( j-^IQi) m + Mg(H))- 9 } 

Remark 7. The gamma exponential tilting operation is quite special and we point out that for 
a fixed value of 9, Yg jC cannot achieve all possible distributions on [0, 1]. So for example when = 1, 
Y\.\ cannot be Uniform[0, 1]. This is equivalent to noting that a density on [0, 1] of the form 

« Jl~y) 

does not exist. In general, the Uniform[0, 1] is possible for Yg iC , only when < < 1. 



3.2 An example with some connections to the occupation time of skew 
Brownian bridge 

Set H = U[0, 1] to denote that the corresponding random variable X is Uniform[0, 1]. It is known 
from Diaconis and Kemperman (1994) that the density of M 1 (C/[0, 1]) is 

(6) - sin{iry)y- y (l - y)~ {1 ~ v) for < y < 1. 

Note furthermore that W — G\Mi(U[Q,l]) is GGC(1, U[0, 1]) and has a rather strange Laplace 
transform, 

E[ e -AGf 1 M 1 (ii[o,i])] =e ( 1 + A )-(^ 1 ). 

We can use this fact combined with the previous results to obtain a new explicit expression for 
the density of what we believe should be an important mean functional and corresponding infinitely 
divisible random variable. 

Proposition 3.3 Let G\/E be the ratio of two independent exponential (1) random variables having 
density ((dx)/dx = (1 + x)~ 2 for x > 0. Now let L denote a GGC{1,Q random variable, with log 
Laplace transform logE[e _AGlMl ^)] = — j^-log(A). Then equivalently the Levy measure of L is 
given by 

Furthermore L = G\M\{Q, where M\(Q has density, 

1 X x 

/mi {x\C) = - sin(7r— — )x for x > 0. 



PROOF. First note that it is straightforward to show that E[e L ] = E[(l + Afi(C)) X ] = e 1 . 
This fact also establishes the existence of L. Now we see that G\/{G\ + E) = f/[o,i]- The result then 
follows by applying statement [(ii)] of Proposition 3.2 to (6)D 

In view of the remarks in Section 3.0.1, we now give a description of the laws of the subordinators 
associated with the two random variables above. 
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Proposition 3.4 Suppose the Z(t) is a subordinator where Z(l) is GGC(1, U[0, 1]), then for < 
t < 1, Z(t) = GiMi([/W[0, 1]), where Mi([/W[0, 1]) has density, 

— sin(7ri(l - y))y t{ ^-^- x {\ - y)~ t{1 ~ v) for 0<y<l. 



IfZ(t) is such that Z(l) is GGC(1,() then then for 0<t<l, Z(t) = GiMi(C (t) ), where Mi(C (t) ) 
has density, 

— sm(n — - — )x<- 1 + ! °)~ 1 for x > 0. 

7T V l+x' 3 



PROOF. This now follows from Theorem 3.1, Proposition 3.3 and (6). □ 

Once we have the density in Mi(() we can then extend the result for Mi(U[0, 1]) to that of 
Mi (Oy 2 p ), where 0^ 2p denotes the distribution of the random variable 



q 2 E' 

In the notation above Bp 1 ^ 2 ' 1 ^ 2 ^ (s) denotes a p-skew Brownian motion. Hence A^ 2? denotes the 
time spent positive by this process up till time 1. The following result is otherwise not obvious. 

Proposition 3.5 Define A br , 2p = p 2 Gi/[p 2 Gi + q 2 E\. The random variable A hr , 2p is equivalent in 

distribution to the time spent positive of a p-skew Brownian bridge having density Oy 2p (dy)/dy — 

p 2 q 2 /[p 2 (l — y) + q 2 y\- Now let L p denote a GGC(1, Oy 2 p ) random variable. Then equivalently the 
Levy measure of L p is given by 

s - 1 E[ e - s/A v^ P ] 

Furthermore L. p = GiMi(Oy 2 p ), where Mi(0^ 2 p ) has density, 

p 2 (i-y) 



irq 2 y 



/m i( o?7, )(y) = — sin y (1 - y) **n-«)+,*v for < y < 1. 

u i / 2 -p' 7r \p {i-y) + qyj 

Where for c = p 2 /q 2 

l/K p = cE[(l + cM{()- 1 }=ce Jo (!+-)=» dx 

Remark 8. The specific densities in Proposition 3.3 and 3.5 yield the not immediately obvious 
identities. 

and, more so, 



f 1 loeClx — v\)v 2 a 2 ( q2y . . p2{1 - v) .. 

$ 2 W = - TYT ^f^rfy = log Kp y->Hi-v) + ,*v (1 - y)~ tHi-vMv 

Jo [p 2 {^-y) + q 2 y] V 

where these are appropriate versions of $. 
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3.3 Reconciling some results of Cifarelli and Mellili 

To further illustrate our point we show how to reconcile two apparently unrelated results given in 
Cifarelli and Mellilli (2000). Let Ai/ 2) i/2 denote the distribution of the arcsine law, that is a Bi/2,1/2 

random variable. Cifarelli and Melilli (2000, p.1394-195) show that for all 6 > M e (A 1/2)1/2 ) = 
^0+1/2,0+1/2- Now define the probability density 

IT 

Cifarelli and Mellili then show that for 9 > 1 M g (g 1 / 2 ) has the density proportional to 

Hjort and Ongaro (2005) recently extend this result for all 9 > and also note the normalizing 
constant appearing in Cifarelli and Mellili (2001) is incorrect. Here however we note that if X has 

density Q\/ 2 then X = Gi/2/G[, 2 where Gi/ 2 and G[, 2 arc independent and identically distributed 
gamma random variables. Now using the known fact that 

d G 1/2 /G' 1/2 

Sl/2 < 1/2 - G 1/2 /G' 1/2 + l' 

we see that £1/2,1/2 is a special case of A\ in Proposition 3.2. It is now evident that one could use 

the result Mg(A 1 / 2il / 2 ) = B 0+1 / 2 0+1 / 2 „ coupled with statement [(ii)] of Proposition 3.2, to easily 
obtain the density of Me(g\/ 2 ) for all 9 > 0. Similar to section 3.2 one could then use the density of 
Mg(gi/ 2 ) to obtain results for mean functionals based on the law of p 2 Gi/ 2 /[p 2 Gi/ 2 + q 2 G' 1 ^ 2 ]. We 
will encounter this class of models again in the next coming sections and see how they arise as laws 
of occupations times of Brownian bridge and related models. 



4 A tilted positive Linnik process 

In this section we present details of a subordinator, Z, which is a FGGC such that Z(l) has the 
Laplace transform, 



(7) 



1 (l + (c + 6A) a ) 



1 + cf 



for parameters < a < 1, > 0, c > and b > 0. The first important thing to note is that Z(t) is 
equivalent in distribution to that of Z{\) with the parameter 9 replaced by 9t. This as we shall show 
more specifically, yields the desirable property of being able to explicitly identify the distribution of 
Z(t) for all t. In this generality we believe that this process has not been studied in any detail, but 
setting c — 0, we see that the Laplace transform becomes 

(i + b a x a y e . 

It is easy to verify that this corresponds to the random variable defined as 

bL a g := bGj S a . 

The random variable, L a fi-, is known in the literature [see for instance Bondesson (1992, p. 38) and 
Devroye (1990, 1996)] and is sometimes called a positive Linnik process and has a host of interesting 
properties and distributional representations. Now let us introduce the random variable which will 

play a key role throughout the remainder of the paper. Let X a = S a /S' a denote the random variable 
having density 

, \ sin(-7ra) 

8 e°(v) - v 2 o^oo — / rn for v > °- 

7r y za + 2y a cos(7ra) + 1 
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Then the Levy exponent associated with L a ^ is 

/>oo 

$ a ,o{w) :=01n(l+u/*) = / (1 - C - Xs )l . a (s)ds 

Jo 

where 

rvft rvG 

hAs) = ^-Ms a ) = ^E[c- sX °] = aBs-^-'l**] 
is the Levy density of the Linnik process. Specifically, 

is the Mittag Leffler function, which equates with the Laplace transform of S~ a . The equivalences 
involving X a are probably not that well known but serve to identify the density of X a (8) as the 
Thorin measure, v of L a j. This Thorin measure is identified by Bondesson ((1992, p. 38) and 
confirms the fact that L Q; g is an FGGC. That is, 

poo />oo 

$ a ,fl(w) := 0a (l-e-^s" 1 / e~ s/x Qa (x)dxds. 
Jo Jo 

One of the things will be looking for are alternative representations of the distribution of L a g. For 
instance it is known from Devroye (1996), that when 6=1 one has L Qi i = GiX a , where G\ is 
exponential (1) and X a has density (8). Summarizing, we see that L a £ is a GGC(a#, g a ) and (7) 
corresponds to the class GGC^ b ' c ^ (a0, g a ) We will establish various distributional equivalences which 
then lead to results for Dirichlet mean functionals and Bessel occupation times. First we describe 
some more pertinent features of X a . 

Proposition 4.1 Let X a = S a /S' a , having density (8). Then, 
(i) The cdfofX a can be represented explicitly as 



(9) F x (x) = 1 - — cot" 1 [ cot(7ra) + 
(ii) Its inverse is given by 

(10) FxAv) = 



sin(7ra) 



sin(7ra(y)) 



_sin(7ra(l — y)) 

(Hi) The equations (9) and (10) yield the quite useful identity, 
(11) sin(7ra(l - F Xa (y))) = y~ a sm{naF Xa (y)) = 



sin(7ra) 



[y 2a + 2y a cos(Tra) + 1] 1/2 

Proof. This derivation of the cdf is influenced by arguments in Fujita and Yor (2006) where it 
becomes clear that it is easier to work with the density of [la] 1 '". The cdf of this quantity is easily 
obtained. Statements [(h)] and [(hi)] then become apparent. □ 

Remark 9. There are several things interesting about the alternative representation of the Levy 
measure involving X a . First, from Chaumont and Yor (2003, sec. 4.19 and 4.21),. the distributional 

equivalence can also be seen to arise from the fact that G^ Q = G\/S' a , where S' a is a stable random 

variable independent of S a , and the fact that X a = S a /S' a . The density of X a can be traced back 
to the work of Lamperti (1958) but arises later elsewhere. One notes also that the density for X a 
has a simple form as compared with a stable law of index < a < 1. 
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Remark 10. One notes that the random variable L a ,g is conditionally a stable random variable 
and hence is heavy-tailed. So this limits its practical applicability to a variety of problems. A natural 
way to create a random variable which has moments is to exponential tilt the density of L a g. This 
is precisely how the Laplace transform, (7), comes about. It is easy to see that, 



E[e- 



-AZ(l)i _ 



E[e" 



-(6A+c)L a , e l 



E[e 



-cL a 



In other words the density of Z(l) is given by 



1 , h c -T(P/9) 1/a 

fz(i)(x) = - b {l + c a )e- xc ' b f La , e {x/b) = b{1 _ p) fr„,.(s/ft). 

4.1 Connection to occupation times of Bessel bridges and P a ,e(C) 

Let Ji > J 2 > • • • denote the ranked jump sizes of an (a)-stable subordinator such that J2^ =1 Jk = 
T a o is a unilateral (a)-stable random variable. Recall from Pitman and Yor (1997a) that the sequence 
of ranked probabilities (P^ — Jt/T a ) is said to have a two-parameter Poisson Dirichlct with 
specification PD(a,0). Furthermore, for a9 > 0, one obtains the PD(a, a9) law of (Pk) by mixing 
over the conditional law of [Pk)\T a fl — t with respect to the distribution of the random variable, say 
T a , a e, which has density f(T a , at g e dt) oc t~ ae F(T a fi e dt). One may then introduce, independent 
of (Pk), a sequence of iid random variables (Zk) having some common non- atomic law. Then as in 
Pitman (1996), see also Ishwaran and James (2001), analogous to the Dirichlct process, one may 
define the class of PD(q, a9) random probability measures as P a , a e{-) = Y^k=i P^Sz k (•)• 

Now it will be shown that the exponential tilting operation to obtain (7) reveals a strong, 
albeit initially unexpected, connection to work of Pitman and Yor (1997b) and Barlow, Pitman and 
Yor (1989), on occupation time models for skew Bessel bridges and more generally laws of P a , a e(C) 
for some set C such that E[P Qi g(C)] = p. 

This is seen for the case where c ^ 0, which allows one to to rewrite (7) as, 

(12) 



(g+p(l + |A) a ) 



for p = c a /l + c a and q = 1 — p. We recognize that for b = c, this equates with the Cauchy-Stieljtes 
transform of order a9 of a two parameter (a, a9) Poisson Dirichlet random probability measure 
evaluated at some set C, denoted as P a , a e{C), where E[P QiQ g(C)] — p. That is, 



E[(l + AP a , afl (C)) 



-0m 



(q+p(l + \) a ) b 



Furthermore, when 9 = 1, this is the Cauchy-Stieltjes transform of order a of the time spent 
positive up to time 1 of a p-skew Bessel bridge of dimension 2 — 2a, as can be seen from Pitman and 
Yor (1997b, eq(75)). This random variable can be represented as, 

< p = Pa AC) - f l{B^ a \s) > 0)ds 
Jo 

where Bp a ' a ^ is the corresponding p-skew Bessel bridge. We also note another important connection 
to Lamperti (1958) and Barlow, Pitman and Yor (1989), Pitman and Yor (1997b). Noting that 
c = (p/q) 1 ^ , the random variable 

(13) ^ p i p„, 0(C) i i i f I(fl«(.) > 0),. 
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equates with the time spent positive by Bp , now a p-skew Bessel process of dimension 2 — 2a 
with skewness parameter p, up to time 1. When p = 1/2, one obtains the usual Bessel processes. 
Lampcrti (1958) shows that the density of A a ^ p is 

(14) A a , p (dx)/dx- Pq Sin(mr) (1 - 



7T [q 2 x 2a + p 2 (l - x) 2a + 2pq x a (1 - x) a cos(a7r) 



We see that setting p = 1/2, a = 1/2, yields Levy's(1939) famous result that the time spent positive 
by Brownian motion up to time 1 has the Arcsine distribution. That is, Ai/2.1/2 is Beta(1/2, 1/2). 
Now interestingly from (9) and (10) we obtain a closed form expression for the cdf and quantile 



function of A atP , as 



(15) F Aa ( y ) = l-— cot" 1 (cot(Tra) + — qy a ) 

p na \ p(l — y) sm(7ra) J 



and its inverse given by 

1/a 

(!6) F t.M = 



q sin(7ra(l— y)) 



Fx Av) 



i + 



p sin(7ra(y)) 



q sin(7ra(l — y)) 



1/a 1 I Tp 1 



Remark 11. Hereafter, for c ^ 0, we shall set c a = p/q. Note that if B(s) denotes in a generic 
sense a Bessel process or bridge, then the interpretation of p is that the p-skewed version of B(s) 
has the property that 

p = P(B(s)>0). 

Now we relate the more general class of PD(a, a6) models to occupation laws of processes and 
their accompanying local times. Noting Pitman and Yor (1992, p. 332) let (4"\i > 0) denote 

the right continuous local time of a Bessel process, and let S a (t) = Sat 1 /" denote an (a)-stable 
subordinator, which satisfies the identities in law, 

< \ ei a) ■ - 1 

(17) S a (s) = mS{t : If' > s} and , -1— 



t<* (S a (s)) a (S a )° 

/ a \ l/a ^ 

For our purposes, we may set S a (l) — T afi , and furthermore [t\ ] = l/T Qj o- Now using the scaling 

property (17) one may construct local times {l\ a ' ae ^ : t > 0} with laws specified by P([£^"' a8 - ) ] ^ e 
ds) — P(l/T Q . Q e € ds). As a special case, £( Q ' Q ) denotes the local time of a Bessel bridge, say 
B^ a \ Associated with {l[ a - ae) : t > 0} are what we shall call , generalized Bessel bridges say 
B (a,afl)( t ) with law specified by P(B^(t)\i[ a) = s)P(e{ a ' ae) E ds). Such processes may be found 
in Definition 3.14 of Perman, Pitman and Yor (1992). Now letting B i p a,a6 " 1 denote ap-skewed version 
of such processes we define their times spent positive up to time 1 as, A^p 00 ^ satisfying, 



A { a a p a9) = P a , aS {C) 



Note that B^ := B^V . 



Remark 12. The representation of the occupation time of a skew Bessel process in (13) was 
given by Barlow, Pitman and Yor (1989). It will play a fundamental role in our understanding 
of the construction of randomized versions and related matters. In effect this boils down to the 
interpretation of c as c a = p/q 
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Remark 13. In this remark we demonstrate how the positive Linnik may be interpreted as the 
distribution of a time changed occupation time. It is known, see Barlow, Pitman and Yor (1989) or 
section 4 of Pitman and Yor (1992), that for S a (t) an inverse local time, 

A a . p (S a (t)) = [ Sa(t) > 0)ds = S a ( P )t^ a = S a (pt) ipV-tVa^ 

Jo 

is an (a)-stable subordinator. Hence letting (Tg(t),t > 0) denote an independent gamma subordi- 
nator satisfying for each fixed t, Tg(t) = Get, it follows that 

A a , p (S a (T e (t))) = S a (pT 9 (t)) = p 1 l a G 1 J t a S a =p l l a S a {T e {t)). 

Remark 14. The random variables T~^e are the a-diversity of the PD(a, a6) exchangeable 
partitions as described in Pitman (2003, Proposition 13). Furthermore it may be read from Perman, 

Pitman and Yor (1992, p. 31) and Pitman (2006) that for 6 > 1 T a ,g a = T a)a e- a Ue a ,i-a, where 

Uoa,i-a is independent of T a $ a but not of T a<a e-a, and marginally Ug aj i- a = B a e t i- a - We will 
sometimes refer to all the relevant local times and occupations times, already defined, as PD(a, aQ) 
processes. 



4.2 Distributional results for the Linnik class and local times 

We first establish various distributional identities related to the positive Linnik random variable. 
Importantly, we will show that random variables based on the PD(a, aQ) models are Dirichlet mean 
functionals. 

Remark 15. Throughout we will be using the fact that if X is a gamma or (a)-stable random 
variable, then the independent random variables X, Y, Z satisfying XY = XZ imply that Y = Z. 
For precise conditions see Chaumont and Yor (2003, sec. 1.12 and 1.13). 

Theorem 4.1 Let L a fi = G X J a S a denote a generalized Linnik random variable i.e. a GGC{a6, g a ). 
Then, we have the distributional equivalences 

(i) For all 6 > 0, L a fi = Gg a S a /T aya g = Gg a M a g(g a ), which implies 

(ii) In particular when 6=1, L a ,\ = G\X a = GiB a _i- a M a (g a ). Hence 

X a = B aA _ a M a (g a ) = MiG?^) 

(in) For < 9 < 1, 

T A T3 l / a n Y A n T3 Set d_ n u S a 

'-'afi — O ei _ e lj\A. a — Lrll30 a; l_0 a — — tsea+l-a&ae^-a—, 

-Lct,a& J-a.aO 

(iv) Statement (i) implies G X J a = G a g jT aM g, for 9 > 0. For < 9 < 1, 

r ,l/a d_ G a g d_ G\Bj s _ g G a e+l-a „ 

^0 — m — o ~~ rp &a0,l-a 

When 6 = 1, then 

fil/a d_ G a G\ 
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Proof. The first equivalence in statement (i) is perhaps the most non-obvious. The result is 
obtained by manipulating the density representation of G^"^ as follows. Let f a denote the density 
of an (a)-stable random variable. Now note that the density of L a fi is obviously expressible as, 

f Lae {y)=Cy ae - 1 / s- 9a e-^ a f a (s)ds = Cy ae - 1 / s - fe E[ e -W s » s "]/ tt ( s )d s . 
Jo Jo 

for some constant C. Now it remains to write 

/>oo 

E[e" (!//s)s °] = / c- vy/s f a (v)dv. 
Jo 

The result is then obtained by algebraic manipulations. The second equivalence is immediate since 
L a fi is GGC(a8, g a ). The remaining statements are straightforward applications of beta-gamma 
calculus and independence. □ 

Remark 16. Statement [(iv)] generalizes the known case of G\^ a = G\/ S a . It is interesting to 
note that this special case can be interpreted in terms of stochastic processes, as, 



Ja) d_ na 1 a n 



where £^ is the local time of a Bessel process with dimension 2 — 2a considered at the independent 
exponential time G\. This description, and further references, may be found in Chaumont and 
Yor (2003, p. 114). It is also known from Barlow, Pitman and Yor (1989) that the local time for the 

bessel bridges, 4"'"'' evaluated at G a , satisfies £^' a ^ = Now using the scaling property in (17) 
it follows from statement [(iv)] that, 

(Go 

(Ge a ) = f = ^ 

See section 5 for an extension of this idea to more general random times. 
Theorem 4.1 establishes the distributional results, 

x a i B aA _ a M a (g a ) i m^) = [(e { £; a ij] 1/a s a i [(4 Q j)] 1/Q 

where (£^) denotes the local time of a Bessel process evaluated at an independent unilateral (te- 
stable random variable, S a . The fact that it is representable as a Dirichlet mean functional, Mi(g^) 
coupled with the explicit density of X a , leads to an important identity below, which plays a funda- 
mental role in obtaining explicit densities throughout the rest of this work. 

Proposition 4.2 Define y a {x) = f™ log(\x - y\)g a (dy) = E[log(|x - X a \)]. Then for < a < 1, 

y a ( x ) = — \ og (x 2a + 2x a cos(aTr) + 1). 
2a 

Note that y a {x) is a special case o/$(x). 

Proof. From Theorem 3.1. it follows that the density of X a = Mi(g^) satisfies the equiva- 
lence, 

Qa ( x ) = l s in(7ra[l - F Xa (x)])e~ s '^ x) x a - 1 . 

TT 

Solving this expressions for 5f a (x), and applying the identity in (11) concludes the result. □ 
With this we obtain explicit expressions for the cdf and density of M a g(g a ) as follows; 
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Theorem 4.2 From Theorem 4.1, M a g(g a ) = [l^'" 6 ^] ^ = S a /T ata g. The form of the cdf for 
M a e(g a ) for all a9 > 0, is given by (1), with 9 := a6 , and 

... 1 sm(ir9aFx (x)) 
A a g(x\g a ) = - 



* [x 2a + 2x a cos(mr) + if /2 

ore, 

a9x a ~ 1 [sin(7ra[l - 6F Xa (a;)]) - x a sm(ir6aF Xa (a;))] 



where Fx a is given in (9). Furthermore, a general expression for the density is obtained from (3) 
with 9 :— a.9 and 

(18) d a g(x\g a ) 



* [x 2a + 2x a cos(aTr) + l] e/2+1 

□ 

Proof. The bulk of the result is a straightforward application of Proposition 4.2 combined 
with the explicit forms of the densities in Section 3. The expression in (18) follows by, noting the 
explicit form of the density g a , and applying the identity sin(w — z) = sm(w)cos(z) — sin(z)cos(w), 
with w — na and z = n9aFx a (x) □ 

From this we obtain the most explicit case of a6 = 1, and the important case where a9 = a, 
related to the local time of a Bessel Bridge evaluated at an independent stable time. 

Corollary 4.1 Consider the random variables in Theorem 4-2 then, for 9 — 1, corresponding to the 

Bessel bridge, the random va 
where d a (x\g a ) simplifies to, 



Bessel bridge, the random variable M a (g a ) = [t^^] ^ = S a /T a , a has density determined by (3) 



, , , , ax a (1 — x 2a ) sin(7ra) 
d a {x\g a ) = — - — r 2a — ; —^2 forx>0 



2x a cos(aTr) + 1] 



When 9 = 1/a, the density of M x {Q a ) = [i^] ° = S a /T atl is given by 

t , i x a / I \ 1 sin(7rFx Q (x)) 

fM 1 {x\g a ) = Ai(ar \Q a ) = -— — -j^. 

K[x 2a + 2x a cos(aTr) + 1] 1/2q 

Proof. Apply (11) to obtain the expression for d a (x\g a ).U 

Remark 17. It is evident from Theorem 3.1 and Proposition 4.2, that we can also obtain similar 
types of expressions, as in Theorem 4.2, for the densities of Mg(g^) = B a g t g(i_ a )M a g(g a ), for all 
6 > 0. For brevity we do not give that here as it can be easily deduced. Instead we concentrate on 
the interesting class of M^g^a^) = B a g t i- a gM a g(g a ) for < a6 < 1. These have the densities of 
non-integral form. 

Proposition 4.3 For < ad < 1, 

M 1 (g^) = B aB ^ a gM a g{g a ) = [{l%£lj] 1/a S a 

and has density 

, 1Q . 1 8m{eC a (x))x ae - 1 

71" [x 2a + 2x a cos(aTr) + lf /2 

where C a (x) — cot -1 (cot(7ra) + sin ^ xa ) ) • The density equates with g a when 6 = 1. For < 9 < 1 
(19) equates with the density of B X J^_ e X a . 
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Proof. This is a consequence of Theorems 3.1 and 4.1 □ 

When a = 1/k, we obtain representations in terms of products of gamma random variables. 

Proposition 4.4 Suppose that a = 1/k for some integer k = 2, 3 . . .. Then since for 9 > 0, Za/fc,0 = 
Si/ k G k g = G g/k M g/k (g 1/k ), one deduces that, 

s k ~ 1 

M e/k (Qi/k) = Tr-^^ = k h S 1/k J] G«+i 
This implies T X j k gj k , is such that 

T J — £kk ii Gs -¥ 

J-i/k,e/k fj[ k 

distribution. One may obtain further representations by using the known factfsee Chaumont and 
Yor (2003, p. 113), that l/S 1/k = k k \\ k Z\ G j/k . 

Proof. The result may be deduced from the equivalence GgSi/ k — Gg/ k M e / k (gi/ k ) established 
in Theorem 4.1 and the following identity 

fe-i 

Gg = (G fc( |)) = k k Ge_ JJ G {g+j)/k , 

which is found in Chaumont and Yor (2003, p. 113). □ 

Remark 18. As mentioned in Section 3, it is known from Cifarelli and Melilli (2000, p. 1394) 

that M s / 2 (gi/2) = G(e+i)/2/Gi/ 2 , for 6 > 2. However the other results in Proposition 4.4 arc 
apparently new for integer values, k > 2, and all > 0. 

4.3 Results for tilted Linnik laws and occupation times 

We now address the case of the tilted processes which connects with the occupation times. 

Proposition 4.5 Suppose that Z is GG& b ' c \a0, g a ) Then, its Levy measure may be read from 
Propositions 3.1. and 3.2 In particular a GGC^ c ^(a8, g) is a GGC(a9,A a ^ p ) with Levy measure 
expressible as 

Additionally for L c a GGC(a9,A a , p ) random variable, we have the following distributional equiva- 
lences 

(i) For all 9>Q, L c = Gg a P a . a g(C) = G Sa M a g(k a , p ), hence M a g(A a . p ) = P a .g(C) 
(ii) As a special case, when 9=1, A b l p = M a (A atP ) 

(m) When 6 = 1, L c = G x A* ap = GiB q4 _ q M q (A q , p ) = GxB a ^ a P a ^{C), where A* ap = 
B aA - a M a (A a . p ) = B a!l _ a P a>ct (C) = GiMi(A^), has density 

(20) ii_^A Q , p (dx) 

Note A* a p equates with the random variable of Pitman and Yor(1997b, Proposition 15). 
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Proof. The correspondence to the P a , a e(G) models follows from (12) as mentioned earlier. All 
other results are consequences of Proposition 3.2. and the use of Beta Gamma calculus □ 

Remark 19. From Pitman and Yor (1997b) and Barlow, Pitman and Yor (1989, p. 307), let 

g a = sup{t : t < 1 : B^ — 0}, where it is known that g a = B a ^_ a . Proposition 15 in Pitman and 

Yor (1997b) and their subsequent discussion establish the results A a ,. p (g a ) = g a A h £ p . with density 
(20). Our result shows this random variable is equivalent in distribution to the Dirichlet mean 

functional, Mi(PSa} p ) = B a ^ a M a (A a , p ). Note that in view of Proposition 3.2. this can be seen to 

arise by first noting again that Mi(g^) = B a! i_ a M a (g a ) = X a where cX a /(cX a + 1) = A a , p and 

then applying the change of measure to get the distribution of M\{A^[j>)- As mentioned in section 3, 
the construction of Theorem 3.1, and indeed Theorem 3.2, were devised in part to understand the 
mechanics of the result of Pitman and Yor (1997b) in more generality. The next result extends this 
idea. 

As mentioned in the previous remark the next result provides an extension of the distributional 
result of Pitman and Yor (1997b, Proposition 15). It also will yield our first concrete example of the 
occupation time of a randomly skewed Bessel process. This may also be seen as a precursor to the 
quantities we will encounter in Section 5. 

Proposition 4.6 ForO < 9 < 1, andGg andGi-g independent, define^ = pBg^_g/([q + pBg^^g}) = 
pGg/(Gg + Gi-gq). Then define 

m) A d cM^) 

a - i; cBl'°_ g X a + l cM.ie^ + l 
The explicit density of the quantity in (21) is given by, 

Ja,e{y) = 



TT [ y 2a q 2 + 2qpy a (l - y) a cos(aTr) + (1 - y) p 
The density of B a g A - a gP a . a g(C) = B a g :1 ^ a gA ( a,p a6) = B a g.i^ a6 M a g{K a . p ) = Mi(A£"p ) ), is 

y ll > 7, -gJa,e{y) jh- 

(! - P) n [y 2a q 2 + 2qpy a (l - yf cos(a^) + (1 - yf a p 2 } 

More generally the results hold for all < Bet < 1, where f a ,e{y) is the density of the random variable 
cM 1 (ei Q0) )/(cM 1 ( (? L a8) ) + 1), and (22) is the density of M^A^P). 



Proof. The density of (21) is obtained from Proposition 4.3. To obtain, (22) note that we are 
using the fact G a $P ata e(C) = GiB a g.i- a gP ata g(C). Now it becomes evident that this distribution 
is obtained from exponentially tilting (p/q) 1 ^ a GiMi(g < ^"°' > ) in the sense of Proposition 3.2 and 
Proposition 4.3. □ 

Remark 20. At this point we could use Theorem 4.2, combined with Proposition 3.2, to obtain 
expressions for the density and cdf of A^j," 6 ^ = P a , a e{G) — M a g(A a ^ p ). Or one could use directly 
(16) and (15). These would provide alternative expressions for P a ,a.6 

(C) obtained in James, Lijoi and 

Priinster (2006). In that work, the authors addressed the case of more general functionals P a , a g(g), 
where they obtained these laws by a direct inversion of the appropriate Cauchy-Stieltjes transform. 
That work does not address functionals such as M a g(g a )- A description of the laws of A^^* 6 ^ will 
appear as a special case of the forthcoming Proposition 5.11. 
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5 GGC/FGGC stable compositions and occupation laws for 
randomly skewed Bessel processes 

Raising things to the level of processes we recall the known fact that a positive Linnik process 
equates to L a ,g{t) = S a (Tg(t)) = G^J^Sa, where in the second equality S a (t) is a stable subordinator 
independent of the gamma subordinator, Tg(t). Based on our previous results we now study the class 

of models S a {Z e {t)) = [Z e (t)] 1/a S a where Zg is a GGC/FGGC. A special case is where Zg(t) is 
itself a stable subordinator of index < f3 < 1, say Sp(i), which satisfies the important identity 

S a (Sp(l)) = [Sfj^^Sa = S a/ 3, an (a/3)-stable random variable. Note that in this case Zg = Sp is 
not a FGGC. 

Remark 21. Obviously compositions of a stable subordinator have been previously studied 
from several important perspectives. For instance, this operation has recently been shown to play 
an interesting role in applications involving coagulation/fragmentation phenomena as described in 
Pitman (2006) and Bcrtoin (2006). 

Here, using the property that Zg(t) = Gg t Mg t (v), we take a different view of S a (Zg(t)) as being 
equivalent in distribution to Ggj Q 5 Q ,(t)[Me t (^)] 1 ^° ! . That is to say, the viewpoint that these are scale 
mixtures of positive Linnik random variables, which leads to a variety of interesting consequences. 
Lets call this class GGC Q (6', v) which will also denote the law of the random variable S a (Zg(\)). 

Now, from the description of the Levy density of a Linnik random variable in Section 4, it is 
evident that the Laplace transform, at t = 1, of S a (Z$(l)), can be expressed as 

E[(l + A a Mj(i/)) - '] = e -*»° (A) 

where 

^ Q) (A) = / 01og(l + \ a r)v{dr) =a6 E[log(l + A^r 1 /")]^*). 
Jo Jo 

That is, the Levy measure of S a {Zg{\)) ± Gl ,a S a [Mg( V )] 1/a can be expressed as 

f>OG 

aOs- 1 / E[c- S ^ x " r )}v{dr) for s > 0. 
Jo 



5.1 Tilting and Randomly Skewed Bessel Bridges 

In parallel to Section 4, we also discuss its tilted version, which as we shall see connects naturally with 
the idea of randomly skewed Bessel processes and bridges. Using the same type of exponential tilting 
as in Section 4 we call the resulting random variables GGCa'°\0, v). Again, provided that c / 0, 

we see that if a random variable L is GGCa' c \9, v), then L = (b/c)L* where L* is GGCa' c \d, v). 
An initial description of the Laplace transform of L* is 

(23) e[(i + (i + a) v>M 9 (^r 9 ] = w 

E[(l + c*Mg(v)y e ] 

where ^"^(A) = ipg a \c(l + A)) — ipg a \c). Now, by setting p a (r) = c a r/l + c a r) and algebra, this is 
equivalent to 

(24) *e Q) (A) = u9 / E[log(l + XA atPa(r) )]u(dr) = a6 / log(l + Xx)K a , v {dx) 

Jo Jo 
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where 

A a .„(dx) = / K a ^ r) (dx)v(dr). 
Jo 

Before examining this quantity further we first formally define, albeit briefly what we mean by 
randomly skewed Bessel processes and bridges and their occupation times. Let £ denote a random 
variable on [0, 1] chosen independently of a Bessel process {B^ (t) : t > 0}. Then the process 

is said to be a ^-randomly skewed Bessel process if |£ = p is a p-skewed Bessel process as defined 

in Barlow, Pitman and Yor (1989). That is, F(B ( *\s) > 0|£) = £. Say that A a ^(t) := f Q I(B^ a) (s) > 
0)ds is the time spent positive of a ^-randomly skewed Bessel process up to time t. Then conditional 
on £ = p, it has distribution A a . p (t) . Equivalcntly, one has, for A a £ := A a ^(l), 

(25) A f = £ 1/a ^a 

Based on (25) we present an interesting special case. 

Proposition 5.1 Define p a := p a /[p a + q a ]. For < /3 < 1, let £ = Ap yPa = c a X /(c a X + 1) 
which is equivalent to £/(l — £) = c a Xp. Then 

Aa,£ — A a p^p. 



Proof. The result follows from the definition in (25) and the easily verified fact that [Xp] x l a X a = 
X a p □ 

^-randomly skewed bridges are defined in an analogous manner where conditionally on £ = p they 
are p-skewed Bessel bridges. Denote the random variable corresponding to the time spent positive 
up till time 1 of such a process as A b J^, which conditional on £ = p, is equivalent in distribution 

to A b ^ p = P a . a {C), having the Cauchy-Stieltjes transform of order a6 in (12). More generally, by 
the usual PD(a, aff) change of measure, we can define A^'^ which conditionally on £ = p equates 

in distribution with A^p"^ = P a ,ae(C). Now we shall describe a sub-class of such processes which 
equates with a random variable described by the transform (23). 

Proposition 5.2 Let A^'Q^ correspond to a ^g-skewed occupation time with having the specific 
density Qg{du\v)/du = k$(1 — u) 6 Qg{du\v)/ du where Kg = 1/E[(1 + c a Mg(i>)) e ] = e^e ( c ) and 

1 u 

Qe(du) = -ifc-Mei-. \v)du 

(1 — u) L — u 

That is Qg is the distribution of the random variable c a Mg(y) j [c a Mg(y) + 1). Equivalently £0 = 
M8(^ c °' c °') and A^'^ satisfies the following properties. 

(i) The Cauchy-Stieltjes transform of order a6 of A^"'^ is expressible as 

[ E[(l + \Atf ] )~ a0 ]Qe(dp\v) = [ (q+p(l + \) a y 9 Qg(dp\v) 
Jo Jo 

and equals (23). 

(ii) Hence, by (24), A<££ 0) = M a6 (A a , v ), and L = G a gA<££ 0) is GGC{a6, A a , v ). 
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Proof. This result is easily verified by an argument similar to the proof of Theorem 3.2. Hence 
we omit the details. □ 



Now define 



sin(7ra) 



(26) Qolv{x) = / ; : Trv(dr). 

y ' * ' K ' 7T J x 2a + 2x a rcos(an) + r 2 v ; 

We can summarize these results in an equivalent manner; 

Proposition 5.3 Let L = [Zg^)] 1 ^ S a denote a GGC a (9,v) random variable then 

(27) L = G ea M a e( Q ){Me(is)} 1/a = G a gM a g(g a . u ). 

where g a>l , is defined in (26). That is GGC a (9,v) is equivalent to GGC(a9, Q a ,„). Now without loss 
of generality if c ^= 0, then set c — 1 and consider the case where L is GGC a (9, i/ 1 ' 1 )) with Laplace 
transform as in (23). Then L = G a gA^'^ = G a gM a g(k a>l/ ) is GGC(a6, A ajI/ ) In general if L is 
GGCfa'^ (9, v) random variable, then its Levy measure can be expressed as, 

/>oc 

a 



poo 

9s- 1 e-^! b ^ s / E[e- s W x « r \v{dr).U 
Jo 



The next result is specialized to the case where v = H. 

Proposition 5.4 Suppose that L is GGC(a9, Qa.n), then there exists a random variable R with 
distribution H, such that R}l a X a has density g a ,H and L is a FGGC with Levy measure 

a9s- 1 E[e-^ Rl/ax ^}. 

Correspondingly L, the exponential tilt of L with c = 1, is GGC(a9, A at n) and has the Levy measure 

a9s- 1 E{e- s / A «^} 

where — £) = R and A a £ = R 1 / a X a /(R 1 / a X a + 1) has distribution A a ,H- Hence in this setting 
^a'if"' ^ M a o(h- a ,H)i and £g = Mg(H ( cQ < c ")) are Dirichlet mean functionals.D 

Remark 22. Bondesson (1992) describes various features of the class which we call GGC Q (9, v). 
See for instance Theorem 3.3.2 of that work which establishes the fact that this class of models are 
indeed GGC. As noted by Bondesson (1992), compositions of GGC random variables with some 
GGC subordinator are not always GGC. An example is the composition of a Gamma subordinator 
with another GGC process. Note that our representation of the Levy measure in terms of X a R}l a 
and our subsequent usage of it appears to be new. 

Remark 23. It is important to note that special cases of such models have already appeared in 
the literature. Proposition 19 of Pitman (1999), in connection with coagulation phenomena, shows, 
with obvious rephrasing, that the law of A b ^ , where £ = £^1-/3, corresponds to a BETA(a — a/3, a/3) 
random variable. Another interesting case where a model having a particular distribution of A*£ ^ 
may be found is in Aldous and Pitman (2004). Specifically, the distribution of a random lengths 1— Tk 
of a T-partition , as described in equation (67) of Aldous and Pitman (2004), may be re-expressed 
as 

F(T fe G dt) = f P(A b a r p G dx)Il k (dp) = P(A^ G dx) 
Jo 
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where Hk (dp)/ dp = (— log(p)) h 1 /(k — 1)! and hence £ = JliLi ^ f° r (^) independent Uniform[0, 1] 
random variables. We now show how to use a construction of Aldous and Pitman (2004) that 
equates, more generally, with the distribution of A^^ 6 \ Extending the definition in Aldous and 
Pitman (2004, p.24), define for < u < 1 

p(a,a8) 

l \ 

Then it follows that T a a C^ = A^C^ and hence, r^" 6 ^ = A^'^ 6 \ 

Remark 24. We note that one could use the results for the laws of A^j?^ = P a , a e(C) in 
James, Lijoi and Prunster (2006), or the form obtainable from Section 4, to represent the laws of 
i^ 9 ' by mixing with respect to a distribution of £. However it is not a trivial matter to simplify 
such expressions. Here we will use a direct approach. 



5.2 Random time changes for Local and Occupation times 

Here we translate the results in the previous section in terms of random time changes of the local 
and occupation times. 

Proposition 5.5 Let (£[ a ' a ^ ;t > 0) denote the local time under the PD(a,a9)law for ad > 0. Let 
Zg(l) denote a GGC(6,v) random variable. Then it follows that 



G ea Ml /a ( v ) 



"""Vw 1 ^ 6)Me{v) = Ze{l) = GeMe[v) 

(ii) A a . p (T a (Zg(l))) =p 1 / Q [(4 Q fl '! ) )] 1/Q 5aM e 1/Q ( I /) ^ P ^ a L a .e. 



PROOF. Statement [(i)] follows from the scaling property in (17), (27), and the fact that 
M a g(g a ) = [£^g' ae ^] • Statement [(ii)] is straightforward. □ 

The result states that we can arrange for the time changed local times to have any GGC distri- 
bution. We present the next result for a simple illustration 

Proposition 5.6 Let (£[ a ' a6 " > ;t > 0) denote the local time under the PD(a,ad)law for a9 > 0. Set 
Mg(v) — Then it follows that 

/ • \ /)(a,o:0) J a, ad) D d „ 

< - r e" £i c.e,e(i-c.) 

(ii) A a , p {T a {G a e))=p 1 l a G)! e a S a 
5.3 Scaling Calculus for GGC a (f?,i^) mean functionals 

We now show that our results establish a rather remarkable type of scaling calculus for mean func- 
tionals. 

Proposition 5.7 IfL is GGCfa' c) (9, v), it follows that GGG i a m (9, v) = GGC(a9, g a . v ) and GGCL M) (0, v) 

GGC(a9 1 K a , v ). That is in the first case L = G a gM a g(g a , v ) and in the second L = G a gM a g(K a ^). 
These points lead to the following results 

(i) M a9 ( Qa , v ) £ M a8 (Q a )[Mg(u)] 1/a 
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(ii) Following Theorem 3.2 and Proposition 3.1, A 1 ^'^^ — M a #(A aiI ,) = Y a g t i, with M = M a g(g a , u ). 

The results apply for arbitrary Thorin measures v but yield results specific to Dirichlet mean Junc- 
tionals by setting v = H.U 

A variation of this result is given below, 
Proposition 5.8 Theorem 3.1 and Proposition 5.7 imply, 
ft) M^j]) = B aeA - a eM a g{ QatH ) i Mi{ e ( £ e) )[M (H)\V a 

(n) M^S) = X a [ Ml (H)] 1/a 

(ii) M 1 {A^)^B ag , 1 _ a eA^f\ 

Notice thatwhenO< 6<l, M^g^J) = X a Bl^_ e [M e {H)] 1/a but Bl^_ e [M e {H)] 1/a = [M^H^f" . 
Hence this is a special case of statement [(ii)].U 

Notice now that 

_ sin(7ra) f°° x°^r 

ea > H W- n J Q x 2 " + 2x<*rcos(an)+r 2 [ } ' 

The cdf of X a R x l a , T a (x) = ¥{X a R}/ a < x), can be represented as 

(28) TJx) = 1 - — [ cot- 1 ( cotina) + ^fi— ) H(dr) 

ira J \ sm(na) J 

which simplifies in the case of a = 1/2 to, 

2 f°° ( fx\ 

^\li{ x ) = — arctan I I H(dr). 

^ Jo V r J 

Proposition 5.9 Let X a = S a / S' a , having density (8). Let R be a random variable with distribution 
H. Then the random variable X a R}l a has density Q a ,H- Define (x) = E[log(|x — X a R}/ a \)I{x ^ 
XR}/ a )]. Then, 

^a(x) = — / log(a; 2Q + 2x a r cos(mr) + r 2 )H{dr) 
2a Jo 

with derivative 

a " [X> J x^ + 2x a rco S (air)+r 2 { >■ 
^1/2^) = Io°° l °S( x + r 2 )H(dr) = J °° log(l + r 2 /x)H(dr) + log(x). □ 

The next result shows that S^^{x) may be seen as the Levy exponent of certain FGGC models. 
This provides another way for further simplification. 

Proposition 5.10 Let R be a random variable with distribution H. Define W a , x = 2Rx~ a cos(7ra) + 
R 2 x~ 2a and denote its distribution as H*. Then, 

R -Ba&*(x) = E [ e -G e/2 M e/2 (H*)^ x -0a 

When a = 1/2 let Hi/ 2 {y) = H(y/y) denote the cdf of R 2 then 

e -e/2y 1 H /2 (x) = E ^ e -± Ge/2 M e/2 (H 1/2 )^ x -e/2_ 

□ 
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5.4 Density formula 

With this we can obtain explicit expressions for the cdf and density of M a e(g a ,H) as follows 

Theorem 5.1 Recall that M a g(g a ,H) = M a g(g a )[M0(H)] 1 ^ a . The form of the cdf for M a e(g a .H) 
for all ctQ > 0, is given by (1), with 9 := a6, and 



A ea (x\g a . H ) = -sin{KeaT a {x))e- ea ^^ 



where T a is given in (28). /S.i(x\g a ^H) is the density of Mi(g a ,n)- Furthermore, a general expression 
for the density is obtained from (3) with := a9 and 

d ae (x\g a . H ) = —e- 9a ^W [ sin (™[l - 0T a (a;)])B*i(aO - sm(irOar a (x))B% 2 (x)] , 



where 



POO q — 1 

BSi(x) = / -5 ^ -, ; ^H(dr) 

a ' lx ' J Q x 2a + 2x a r cos(cra-) + r 2 K ' 



and 

„2a-l 



roc X 

B" 2 (x) = / — : r ^H(dr). 

' Jo x 2a + 2x a r cos(aTr) + r 2 v ; 



□ 

Now applying Theorem 3.1 we obtain, 

Theorem 5.2 For < a6 < 1, the density of M^g^J) = B aeA _ ae M ae (g a )[M e (H)] 1/a is given 
by 

- sin(7rfe[l - T^e-^^x^- 1 
When 6 = 1, the density of M x {q^ h ) = X a [Mi(H)] 1 ^ a can also be expressed as, 

sin(Tra) [°° x a ~ x r 
^^ = —^1 x 2 « + 2x«rcos(an) + r 2 fM ^ H)dr - 

□ 

Then next results are for the ^-skewed occupation times. 

Proposition 5.11 Assume without loss of generality that c = 1. Then, for all a6 > 0, the density 
of (A Q ,£f) = A^'^ can be expressed as 

- y) ae ~ 2 fM a9 (Y^j\e a ,H) 

PROOF. This is an application of Proposition 3.2. See also statement [(h)] of Proposition 5.7. □ 

The next result is a generalization of Proposition 4.6. 

Proposition 5.12 Assume without loss of generality that c = 1. Then, for all < a9 < 1, the 

density of M\{PS^) = B a 0^^ a gA^£ ^ can be expressed as 



ir 1 - y 
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PROOF. Apply Proposition 3.2 to the density of Mi{g^ 6 ^) given in Theorem 5.2 □ 

The next result is also immediate from Theorem 5.2 
Proposition 5.13 Set £/(l — £) = M\(H), then the density of 

[Mi{H)f a X a 



[M 1 (H)] 1/a X a + l 
is 

- sin( TO [l - T a (-^-)])e-^^) y a-l (1 _ ya-l 

7T 1 - y 

□ 

Remark 25. We mention briefly the following example, in the case of a = 1/2. Choose R 2 = 
Gi/E, then R 2 X 1 / 2 has the distribution of a Pareto TYPE III law with cdf 

and e ~ e/2 ^" 2 ^ = E[e~^ Ge / 2Me / 2 ^]x~ e / 2 = i^i Where the last expression is obtained from 
Proposition 5.10 and Proposition 3.3. 

5.5 Identities for PD local times and skew Bessel Bridges 

We close by showing how the scaling properties discussed in Proposition 5.7 and 5.8, coupled with 
the identity X a p = X a [Xp] l / a , translate into interesting identities for skew Bessel bridges and 
corresponding local times. The first result is an immediate consequence of Propositions 4.2 and 5.7. 

Proposition 5.14 For < (3 < 1, set R = Xp, that is H — Qp, then for 6 > 0, 

M a6 (g a p) i M a e(g a )[M e (gp)} 1/a . 

Equivalently, 

(i) / Q/3 < Qe ) £ p(0,0)rp(c<x8)f 

(it) In terms of the (a) -diversity, 



S a p S a 



Sp 



L P,6 



1/a 



T a p,a8 T a , a 

which implies that T a p^ a g = Ta^elTp^} 1 ^ .O 
Now applying Proposition 5.8 leads to, 

Proposition 5.15 For < (3 < 1, set R = Xp, that is H = gp, which corresponds to g a ,H = g a p, 
then for < a.6 < 1, 

(29) M 1 ( e ( 3 ) )±M 1 (eW)[M e (e fl )] 1/a 

Recall that [£^°^] ^ = X a = Mi(gif') = B a .i- a M a (g a ) , then as special cases of (29), with the 
choice of 6 = [3 and 6 — 1 respectively, 
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(i) X aP = Mi(g^p) = M^Q^MpiQp)] 1 ' 

(ii) M 1 {Q<$)±X a [M 1 {Q fi )] 1 ' a .n 
The next result shows how to recover p-skew Bessel bridges from ^-skewed bridges. 

Proposition 5.16 For < ft < 1, set R = c a Xp, where c — (p/q) 1 ^" , following Proposition 5.2 

this means that A aH = A a p. p . Now define p a = p a / (p a + q a ). Then, £ e = A^ p °^ = Mg(Ap^ Pa ) and 
for this choice, 

(30) I A<$£*> ± M a g(A af3 , p ) = P a(i , a e{C) 

Results recovering p-skew Bessel bridges follows as special cases of (30), with the choice of 8 = ft 
and 9 = 1 respectively, 

fx ) AcaP) A A {*P^fi) 
' v a,l/3 ~~ a/3,p 

(ii) 4( Q ' Q ) — A {af3 ' a) 

Proof. We present a proof for clarity. First, as noted, the choice of R = c a X{j coupled with 
proposition 5.2 shows that A a n = A a p tP , leading to (30). What remains is to verify that £g = A^^\ 
Noting Proposition 5.4, & = Mg{H^ a ^) , but here ij( cCV > c °) = A^, when H = g f3 . □ 

So in closing we see, as a special case, that occupation time of a Brownian bridge up to time 1, 
say Ay 2 := Ay^i/2 ^ ^ ^1/2,1/2 is equivalent in distribution to the time time spent positive up to 
time 1 of a A^f^ randomly skewed process B^ 1/2 \t) for any a = 1/(2 ft) and 1/2 < ft < 1. That 
is to say, this process is randomly skewed by the random variable corresponding to the time spent 
positive up till time 1 of a Bessel bridge of dimension 2 — 2/3 < 1, _B^f ,/3 '(t). 
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